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Abstract 

We consider the massive scalar-tensor theory in the Jordan frame = K^^^ 

and U{^) = (l/2)m^$^, where -F($) corresponds to a constant Brans-Dicke parame- 
ter cobd = l/4i^^. The constraint of the Solar System experiments is < (1/400)^. 
For dustlike matter in a spatially flat homogeneous isotropic universe, we reduce the 
equations of motion to a system of two differential equations of first order which can 
be exactly solved. We obtain simple and explicit expressions for and -^j^ that 
depend only on two parameters, and ^mfl- For K < 1/400 the expansion rate 
H{z) can be practically superposed on the ACDM solution H\{z), up to high redshift 
z, but the equation of state wde{z) of the dark energy is not constant: it presents a 
very slight crossing of the phantom divide line w = —1 m. the neighborhood of z = 
and becomes very slightly positive at high redshifts. 

PACS: 98.80.Cq, 04.20.Jb, 98.80.Es 



1 Introduction 

Supernova observations [1], [2], [3], jl] provide strong evidence for a late-time accelerated 
expansion of our Universe. These data, combined with the observations of the microwave 
background and gravitational clustering, strongly suggest that our universe is spatially flat 
and composed of about two-thirds of an unknown form of negative pressure unclustered 
matter called Dark Energy (DE) |5], [6], [7]. The Friedmann- Robertson- Walker equations 
with a cosmological constant A as source of dark energy are the simplest model (called 
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AC DM), which fits remarkably well with the current cosmological observations [8]. How- 
ever it is difficult to explain the magnitude of A, which is unnaturally small (e.g. [9]). 
In the ACDM model the energy density of the dark energy pde = and the pression 
Pde = —PDE are constants, so that the value of wde = is exactly —1, that is wde is 
exactly on the phantom divide line. The uncertainty on the observations tells us only that 
wde{z) is close to —1 and can eventually cross the phantom divide line at a late time. So 
there is intense activity on models likely to cross the divide line; among them, the class of 
the scalar-tensor theories, which is the subject of this paper. 

Since there is no obvious reason to choose a particular model, defined by the func- 
tions and in many works on the subject the phenomenological aspects are 
emphasized, i.e. on reconstructing these functions from observations. The knowledge of 
the luminosity distance Di{z) and the linear density perturbation 6m{z) is sufficient to 
reconstruct the full theory [10] . It is also possible to consider a restricted class, by giving 
an additional constraint a priori, for example by giving and reconstructing 
from the knowledge of the luminosity distance {Di{z) [TT] or by giving a parametrization 
of W£)e{z) or H{z) and reconstructing and [12], 

In this paper we consider the model = K'^^'^ and f/($) = (l/2)m^$^. Its 

particular appeal lies in that the Brans-Dicke parameter ubd = / (^)^ has a constant 
value 1/4K^ and f/($) is the standard mass potential. From the Solar System experiments, 
i^BD > 40000, hence < (1/400)^ is a small parameter. For dustlike matter in a spatially 
fiat, homogeneous and isotropic universe, we obtain an exact solution for the expansion 
rate H{z), compatible with its ACDM counterpart up to a large redshift. 

Let us note that this particular scalar-tensor model has been studied several times in 
different contexts. In particular, under the same conditions, that is for dustlike matter in 
a spatially flat homogeneous isotropic universe, K. Uehara and C. W. Kim |14] have solved 
this system (in the Brans-Dicke form) and given an exact solution for the scale factor a{t). 
But this function of time is not directly adapted for comparison with experiments which 
require functions of the redshift z, and it seems very difficult to extract an exact H{z) from 
their solution a(t). 

In section 2, we introduce the equations of motion for a scalar-tensor theory in a Jordan 
frame for dustlike matter in a spatially flat, homogeneous and isotropic universe and carry 
out some preliminary calculations that set the stage for our approach. In section 3 we 
derive, from the equations of motion, a system of two differential equations of flrst order 
for u{z) = {H{z)/Hq)'^ and S{z) = ^^^^ where the independent variable is the redshift 
z . In section 4, we obtain an exact solution of this differential system. The quantities 
h{z) = and are explicitly expressed as functions of the two parameters K'^ and 
Qm,o and h{z) is compared with its ACDM counterpart. In section 5, we look at the 
equation of state wde{z) of the dark energy, which is not constant and we show that it 
displays a slight crossing of the phantom divide line w = —1 in the neigbourhood of z = 
and tends to ^K"^ + 0{K'^) at inflnity. We show that the dark energy is negligible at high 
redshift. In section 6, by a local dynamical analysis, we show that the orbit of this solution 
approaches a stable proper node. In section 7, we conclude and explain why it would be 
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interesting to fit the solution H{z) to the data of the supernovae observations. 



2 Scalar- Tensor Gravity and ACDM 

Let us consider a universe where gravity is described by a scalar-tensor theory. The La- 
grangian density in Jordan frame is given by 

L = i {F{<I>)R - (?^^9/i$az/<l> - 2Um + LUg,u), (1) 

where L^ig^u) describes the dustlike matter. 

For a spatially flat homogeneous isotropic universe whose metric is given by 

ds^ = -dt^ + a^{t)dx'^, (2) 

the equations of the motion are 

?>FH^ = Pm + - ?>HF + U, (3) 

- 2FH = pm + ^^ + F-HF, (4) 

dF , • d\J , , 

$ + 3i/$ = 3— Off + 2if2 - 5 

where H = ^ describes the evolution of the expanding universe and pm is the energy density 
of the dust matter. It is useful to have also the equation of conservation of matter 

Pm + SHpm = 0, (6) 

which is not independent of the other equations. 

It is convenient to introduce the dimensionless quantities 

h = ^ , t' = Hot (7) 

where Hq = H{to) is the Hubble constant. The quantities evaluated today at to will be 
denoted with a zero index. 

We successively ehminate pm from the equations ([3]) and (jlj), and $ from the result- 
ing equation and the equation (j5]), then writing these equations, using the dimensionless 
quantities ([7]), we obtain a system of two differential equations for h{t') and 

hHW^-yL)^-h'(l+J-^) + -(-+F.,. ) = 0, (8) 



3 



$" + 3/i<l>' + ^-3F$(/i' + 2/i2) = 0, (9) 



'0 



where $' = ^ etc... Another equation will be useful, from at t'^: 



where 



As indicated in the introduction, we choose now 

The equations ([8]), ([9]), ffTOj) become 

(1 + + + + (gij + ?)|; - ^/^ - i(4,„S + ^) = 0, (13) 

$" + 3/1$' + - 6K'^^{h' + 2/1^) = 0, (14) 

with 

m 

mo = — . (16) 

-"0 

It will be interesting to compare these equations to those of the corresponding ACDM 
model that we recall here: 

= + -, (17) 

a SnGpm , A 

(18) 



a 6 3 



From ffTTl) and ffTSl). we obtain 



+ 2^/ _ ^^^^ ^ 0, (19) 
1 = Qm,o + f^A.o, (20) 



with 



^W-^^ ' ^^^'°-3^- ^^^^ 

Let us note that the convention chosen in the Lagrangian ([T]) implies 8ttG = 1/Fq to 
very high accuracy. 
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3 Reduction to a system of two differential equations 
of the first order 



The redshift z is given by 
From f l22p and ([7j) we have 



1 + ^ = ^. (22) 
a 



= -n ^TTTT' (23) 

{l + z)h{zy 



and the differential equations (1131) and (fT^ are rewritten in terms of the variable z: 



h\l + m') - -(1 + 6ir^)(l + z)h— + (— - + -)(! + zYh 



+3(1 + ^)^^^ - q(^< + = 0' 

'1+^) (''dT rfj + ''^faT ) -2(l + --)ft^^+™S*+6A-* ((1 + .)/.^ - 2/,^j = 0. (25) 
We now express the system (12^ . (125|) as a system of the first order in the new functions 

u{z) = h\z) , ^(.) = ^: (26) 
-1(1 + 6K^)(1 + + (1 + + ^(1 + + + + ^)'^' «W 



;i + z)^5(^)+3ir^(l + z) ) ^+ 



(27) 
(28) 



-12i^2 - 2(1 + z)S{z) + (1 + zf + S\z)^ ^ u{z) +ml = 0. 

4 Exact solution and comparison with ACDM 

We remark that the two equations (12 7p . (12 8 p of the differential system are written in the 
same form as first order linear equations in u{z): 

Ai— + Biu{z) + d = , 1 = 1,2, (29) 



where the coefficients Ai,Bi,Ci depend in general on z,S{z), and on the parameters 
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In order to have a solution u{z) of this system it is sufficient to find S{z) which ensures 
the proportionahty between the coefficients: 



that is: 



A2 = /xA , B2 = fiBi , C2 = /iCi, (30) 



-(1 + z^Siz) + 3K\1 + z))= + 6K^)il + z)) (31) 



-12K^ - 2(1 + z)Siz) + (1 + zf + S\z] 
^^[l + 8K' + \{l + z)S{z) + + + ^?S\z) 



(32) 



"^^ = M^)(4 + ^) (33) 

Fortunately, the system of equations (13T!) . (!32|) and (!33l) has a solution (/x, 5'(2;)). 
The solution of (l33l) gives 

that we insert in (1311) . which is an algebraic equation in S{z). We obtain 

Inserting flMl) and fl35l) in the differential equation (|32|) . we find that it is identically verified. 

We remark that the solution fl35l) of ^(z) is independent of mo and that there is not 
any unknown constant of integration since it is obtained by solving an algebraic equation. 
We have 

= -(TTiin- (36) 

Now, we can obtain the solution u{z) by replacing S{z) given by (l35ll in one of the two 
equations of the system fl27|) . (1281) . We obtain a linear first order equation in u{z): 

fill Tfi ^ 

(1 + AK^){1 + QK^){1 + z)—-{?> + UK^ + 96i^^)M(2) + ^(1 + 47^^)2 = q. (37) 

It will be better for comparison with ACDM to introduce the parameter Qrn,o instead of 
rriQ. We eliminate ttIq from f|T5|) 

mg = 6is:^ [1-^^,0-25(0)] -5^(0) (38) 

where we have used 

|; . _(^?^ . ^5(0). (39) 
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and with ( 136|) we write (137|) in the form 

nil 

(l+4ir2)(l+6ir2)(i+^)^_(3+34ir2+96ir^)M(2)+3+34ir2+96i^^-3(^„,o(l+4K')' = 0, 

(40) 

whose solution is determined by the initial condition m(0) = 1 : 

«(z) = l-fio + ^^o(l + ^) (41) 

where 

n _n 3 + 24K2^48ir4 

3 + 347^2 + 967^4- ^^^^ 

In the ACDM model, equation (fT9|) becomes, with the use of (l20l) . 

- 1(1 + z)^ + «(^) - (1 - ^]^,o) = 0, (43) 
whose well-known solution is 

Ma(^) = 1 -^^m,0 + ^^m,o(l + ^)^- (44) 

Now we know, as indicated in the introduction, that we must have < 1/(400)^ from 
solar experiments. So, it is obvious from (HT]) . and that solution u{z) is very close 
to ma(z). 

From the exact solution fHTl) . fl35|) and fl26|) we have immediately 



/i(^) = V 1-^0 + ^0(1 + ^)^^ (45) 



Ho 

and 

[1 + z) . (46) 



$0 

We illustrate in Figured! for = 1/(400)^ and Qm,o = 0.3, the proximity of the expansion 
rate H{z) with H\{z) the expansion rate of the ACDM model by plotting the relative 
difference {H{z) — Hj\,{z))/ H/^{z) which is independent of Hq. 

Let us observe from ( !38l) and ( !36l) that we have an algebraic relation between mg, K"^ 
and flrn,o which can be written 

ml = QK^{l-nrn,o) + 0{K^). 

This relation shows that the mass m = mo-ffo is very small, which imposes a strong 
constraint on the model. 

It is interesting to give an approximation of the exact solution u{z) to first order in K^: 

UapproA^) = l-^^„,o(l-yi^') + ^^™,o(l-yi^')(l + ^)' + 4i^'^^„.,o(l + ^)'ln(l + z). (47) 
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We see that we can neglect = 1/(400)^ in the two first terms of ( H7|) but the last term 
can dominate for very large z. However this term is still 1/100 of the last but one, that is 
AK^ ln(l + z) = Y5o5 ^ — 5-22 x 10^^^ which is a very large redshift. 

For z = 3000, which marks approximately the transition between the nonrelativis- 
tic matter and the relativistic matter (not taken into account in this model), we have 
AK"^ ln(l + 3000) ~ 2 x 10~^; the last term is negligible and we can write 

UAiz) ~ Uapproxiz) ^ u{z). (48) 

In the next section, we shall return to this approximation, in relation with the equation 
of state parameter of the dark energy wi:>e{z) ■ 



5 Equation of state parameter of the dark energy 

We can define the DE density pde and pression Pde by writing the equations ([3]) (j4]) in 
the form [15] 

SFoH'' = pm + PDE, (49) 

- 2FoH = pm + PDE + Pde- (50) 
With these definitions the usual conservation equation applies: 

Pde = -^H{pde + Pde)- (51) 
And if we define whe = we obtain 

PDE ' 

_ i{i + zy-^-u{z) 

- u{z)-n^,o{i + zr- ^''^ 

Details are given in [15]. From (152|) and (jHj), we write 

, , -1+^0 + 1^^0(1 + .)^+!^ 

wde{z) = , (53) 

i-no + fio(i + zf^T^ - ^]„,o(l + z) 



\3 



hence ^ 

We can also see that we have a slight crossing of the phantom line in the neighborhood of 
z = 0. The condition weie{0) < —1 is written 

3(^)0 -^^m,o<0, (55) 
which is easily verified using Uapprox{z) given by fj47j) : 

- 2K^nm,o < 0. (56) 
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The plotted Figures ([2], [3]) show wde{z) at different scales. Figure (jl]) displays the evolution 
of Wr)E{z), when K decreases from K = 1/400, Woe{z) takes more time close to the value 
-1. 

It is now customary to write u{z) [15J as a sum of terms of nonrelativistic energy and 
dark energy: 

u{z) = fi„,o(l + zf + (1 - n^,o)e{z), (57) 
where e{z) is expressed as a function of wde{z) by 



e[z) = exp 



^^A + woe{z') 



{51 



1 + z' 

In order to evaluate e{z), it is simpler to compare f lFTI) with f H7|) . which is written: 

Uapprox 

(2;) ~ 1 - ^Imfi + ^^m,o(l + -2)^ + 4K^i7m,o(l + z)^ ln(l + -z). (59) 

(It is legitimate to neglect = (1/400)^ in the two first terms.) 
Hence, a good approximation of e{z) is given by 

e(z)-l + 4ir^ ^"!;° {l + z)^\n{l + z). (60) 

It is immediate to incorporate the result of the end of the last section: The part of 
the dark energy (1 — Qrn,o)^ becomes preponderant for very large z but it is yet 1/100 of 
the term of nonrelativistic matter fim,o(l + z)^ for the huge redshift z ^ 5.22 x 10^^'^ . 
At z = 3000 it is perfectly negligible. So this term should be significant only at the very 
beginning of the Universe. But we have not included in our model the relativistic matter 
~ (1 + z)'^ which becomes significant at about z = 3000 and can overtake the dark energy 
~ (l + 2;)^ln(l + z) at large z. This last remark suggests that the dark energy of this model 
is never significant at large redshifts. 

The deceleration parameter q = — ^ is given by 

q{z) = -l + -{l + z)-^y (61) 

The plot of —q{z) (Figure [5]) shows that the acceleration of the Universe in this model 
{flmfi = .3,K = 1/400) begins at z ~ 0.65. It can be verified that q{z) is virtually 
independent of K where K is compatible with the Solar System experiments. 



6 Dynamical study around the exact solution 



The system of differential equations (1271) , fl28|) is transformed to an autonomous dynamical 
system by setting 

v = S{z){l + z), (62) 
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and the change of variable 

[p is growing as a function of time.) 
We obtain: 



P 



\n{l + z). 



(63) 



du + AK^)ml - u{p) {6K^ + A8K^ + AK%{p) + (1 + AK^)v\p)) 

dp ~ 2K^{1 + QK^) ' 

dv _ {2K^ + v{p) + AK hjp)) {ml + u{p) {QK^ - 12K^v{p) - v\p))) 
dp 



4^2(1 + QK^)u{p) 



The exact solution is written 



v{p) = 



1 + 



u{j>) = 1 — ^0 + ^0 exp — 

Let us note that p = corresponds to z = 0. 
There are two hyperbolic fixed points 



(3 + 



)p 



(64) 
(65) 

(66) 
(67) 



Ui 



Ami 



3(3 + 167^2) 



U2 



2\2 



mn 



^1 = 



V2 



-2K' 



(6J 



(69) 



2K^{3 + 3AK'^ + 96K^) ' 1 + AK'^ J ' 

The eigenvalues of the associated linear dynamical system allow the classification of these 
fixed points. 

The first one, (mi, f i), is a saddle point (two eigenvalues of opposite sign) which does not 
concern our solution. The second, {u2,V2), is a stable proper node (two identical negative 
eigenvalues), and the orbit of the exact solution fl67t [661) is the horizontal line of the phase 
portrait which approaches this node. 

The stable node ensures the stability of the exact solution. Let us illustrate this result. 



By taking initial conditions on f (0) in the neighborhood of V2 



-2K' 



, (always u{0) 



by definition), we have given numerically the local phase portrait around the horizontal 
orbit of the solution (1671 166|) (cf figure |6]) . 



7 Conclusions 

We have obtained an exact and simple solution H{z),^{z) for an attractive scalar-tensor 
model, considered several times previously in the literature. If the constraints of the Solar 
System observations are respected {K"^ < (1/400)^), H{z) cannot be distinguished from 
H\{z) up to a very large z for the same fim,o- But the equation of state wde{z) is different 
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from the scaling solution Wde = —1 of the AC DM model since it presents a very slight 
crossing of the phantom divide line w = —1 in the neighborhood of 2; = and tends toward 
(4/3) /s:^ at infinity. 

The dependence of the solution H/Hq in Eq. ( I45l) on the two parameters Qm,o and 
could be fitted with the supernovae experimental data using the maximum likelihood 
method, which involves the minimization of the function 



i=i ^ 



where N is the number of the observed SNe luminosity distances -DL(-2j)obs) the corre- 
sponding errors and Dl{z) = {1 + z)/Hq dz' /h{z') is the theoretical luminosity distance. 
This fit would be at least as good as the fit of the ACDM model which minimizes x^i^mfl, 0) 
since the ACDM solution is obtained for K"^ — > 0. It would be interesting to verify if this 
fit respects the solar condition < (1/400)^. 

Finally, since the history of the expansion of this model is virtually identical to that of 
ACDM, let us see if the observation of the growth of the perturbations could distinguish 
between these two models. The equation of the fractional density perturbation 6m{z) of 
this model pi)j and of its ACDM counterpart can be written 



Mz)Ciz) + - ^) S'^iz) ^ ^(1 + z)n^,oSUz). (71) 

The coefficients of these linear differential equations are known and numerically close since 

u{z) ~ u/^{z) and = (1 + 2;) i+^f^- i until a large redshift. So it seems that it would 
be difficult to discriminate between these two models on the basis of these perturbation 
equations. 
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Figure 2: wde{z)- K = 1/400, ^mfl = 0.3. Slight crossing of the phantom divide. 
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Figure 3: wde{z). K = 1/400, Qmfi = 0.3. 
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Figure 5: acceleration —q{z). K = 1/400, ^2^,0 = 0.3. The acceleration begins at z ~ 0.65. 
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Figure 6: Local phase portrait around the horizontal orbit of the exact solution (dashed 
lines). The continuous lines are the phase portrait of the linearised system. K = 1/400, 
^m,0 = 0.3. 
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